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Abstract 

A finite rank separable approximation for the particle-hole RPA calculations 
with Skyrme interactions is extended to take into account the pairing. As 
an illustration of the method energies and transition probabilities for the 
quadrupole and octupole excitations in some O, Ar, Sn and Pb isotopes are 
calculated. The values obtained within our approach are very close to those 
that were calculated within QRPA with the full Skyrme interaction. They 
are in reasonable agreement with experimental data. 



I. INTRODUCTION 



Many properties of the nuclear states can be described within the random phase approxi- 
mation (RPA) . Among many microscopic nuclear models aiming at a description of the 
properties of nuclear excitations the most consistent model employs an effective interaction 
which can describe, throughout the periodic table, the ground states in the framework of 
the Hartree-Fock (HF) approximation and the excited states in the small amplitude limit of 
the time-dependent HF, or the random phase approximation (RPA). The Gogny's interac- 
tion M and the Skyrme-type interactions are very popular now. Such models are quite 
successful not only for predicting accurately nuclear ground state properties f7|,|8|] but also 
for calculating the main features of giant resonances in closed-shell nuclei PJlOl and single- 



particle strengths near closed shells |TT| . Taking into account the pairing effects enables one 
to reproduce also many properties of collective states in open-shell nuclei [|E^-[l5j . 

It is well known that due to the anharmonicity of vibrations there is a coupling between 
one-phonon and more complex states HQ. The main difficulty is that the complexity of 
calculations beyond standard RPA (e.g., for studying damping mechanisms of collective ex- 
citations) increases rapidly with the size of the configuration space and one has to work 
within limited spaces. From another point of view more phenomenological models that 
assume some simple separable form for the residual nucleon-nucleon interaction while the 
mean field is modelized by an empirical potential well allow one to calculate nuclear excita- 
tions in very large configuration spaces since there is no need to diagonalize matrices whose 
dimensions grow with the size of configuration space. The well-known quasiparticle-phonon 
model (QPM) of Soloviev et al. |4j] belongs to such a model. Very detailed predictions can 



be made by QPM for nuclei away from closed shells [[T(| . 

The possibility to solve easily the RPA problem in a large configuration space when 
the residual particle-hole (p-h) interaction is separable was the motivation for proposing in 
our previous work |T7j a finite rank approximation for the p-h interaction resulting from 
Skyrme-type forces. Thus, the self-consistent mean field can be calculated in the standard 
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way with the original Skyrme interaction whereas the RPA solutions would be obtained with 
the finite rank approximation to the p-h matrix elements. It was found that the finite rank 
approximation reproduces reasonably well the dipole and quadrupole strength distributions 
in Ar isotopes. 

In the present work, we extend the finite rank RPA calculations to take into account 
pairing effects. As an application we present results of calculations for low-lying 2 + and 3~ 
states in some O, Ar, Sn and Pb isotopes. This paper is organized as follows: in Section II 
we sketch our method for constructing a finite rank interaction for the quasiparticle RPA 
(QRPA) case. In Section III we discuss details of calculations and show how this approach 
can be applied to treat different multipole states in wide excitation energy regions. Results 
of calculations for characteristics of the quadrupole and octupole states in some nuclei are 
given in Section IV. Conclusions are drawn in Section V. 



II. HAMILTON! AN OF THE MODEL AND QRPA 

We start from the effective Skyrme interaction p| and use the notation of Ref. [TE 



containing explicit density dependence and all spin-exchange terms rather than the original 
form of Ref || where density dependence at the HF level was introduced by a three-body 
contact force and where some spin-exchange terms were dropped. The exact p-h residual 
interaction V res corresponding to the Skyrme force and including both direct and exchange 
terms can be obtained as the second derivative of the energy density functional with respect 
to the density flnj. Following our previous paper [17]] we simplify V res by approximating it 



by its Landau-Migdal form in the momentum space: 



V res (ki, k 2 ) = Aq' 1 £ [Fi + G m ■ a 2 + (F/ + G\ai ■ a 2 )n • r 2 ] Pi(^r^), (1) 

1=0 K F 

where kj, Oi and Tj are the nucleon momentum, spin and isospin operators, and 
Nq = 2kpm* /TT 2 h 2 with \zf and m* standing for the Fermi momentum and nucleon effective 
mass. For Skyrme interactions all Landau parameters with I > 1 are zero. Here, we keep 
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only the I = terms in V res and in the coordinate representation one can write it in the 
following form: 

Vresin, r 2 ) = N" 1 [F Q (n) + G G {r x )a x ■ a 2 + (iftn) + G^n)^ • a 2 )n ■ r 2 ] 6{n - r 2 ) (2) 
The expressions for Fq,Gq,Fq, G' in terms of the Skyrme force parameters can be found 



in Ref. []n|. Because of the density dependence of the interaction the Landau parameters 
of Eq.(0) are functions of the coordinate r. In what follows we use the second quantized 
representation and V res can be written as: 

Vres = \Y1 ^ 12 34 : dt a 2 a i a 2, ■ (3) 
^ 1234 

where af (ai) is the particle creation (annihilation) operator and 1 denotes the quantum 
numbers (nilxjirrii), 

Vi 234 = y0t(r 1 )0;(r 2 )Ke S (ri,r 2 )03(r 1 )04(r2)dr 1 dr 2 , (4) 



V1234 = J 2 (-) K (ji m ih -m 3 \ J - M)(j 2 m 2 j 4 - m 4 | JM)V^ 2U , (5) 

JM 

where K = J + j$ + j± — M — — and 

^1234 = 0*1 1 1 Y J I lis) (j2 \\Yj\ \U) I M (JU2J3J4) - 

(h\\TjL\\h){j2\\T JL \\j4)h{jij2hji)- (6) 

L=J,J±1 

In the above equation, lj'3) is the reduced matrix element of the spherical harmonics 

F Jm , j = V2J + 1, T$l(r,a) = [Y L x a]f and luijxhhji), hijihhk) are the radial 
integrals: 

r°° dv 
iMihhhh) = J o (Mr) + Fo(r)ri ■ r 2 ) «i 1 (r)M i2 (r)u i3 (r)u i4 (r) — , (7) 

roc (j^if 

Isihhjzk) = J o ( G o(r) + G'o^ri ■ r 2 ) u h {r)u h {r)u h {r)u h {r) — , (8) 
where the radial wave functions u(r) are related to the HF single-particle wave functions: 
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n 



(9) 



As it is shown in [[L7[] the radial integrals can be calculated accurately by choosing a 
large enough cutoff radius R and using a iV-point integration Gauss formula with abscissas 
and weights r k , w k . Thus, the residual interaction can be presented as a sum of N separable 
terms. 

So we employ the hamiltonian including an average nuclear HF field, pairing interactions, 
the isoscalar and isovector particle-hole (p-h) residual forces in a finite rank N: 



H = E I E ( E j ~ K)a' jm a jm 



1 N 

k=l q=±l Xfi 



K (M, k ) + gK (M, k) 



< )+ (r)M^(gr) : + 



E 

L=A,A±1 



(S,fe) . 



(10) 



We sum over the proton(p) and neutron(n) indexes and the notation {r = (n,p)} is used. 
A change r <-» — r means a change p ^ n. The single-particle states are specified by the 
quantum numbers (jm), £j are the single-particle energies, A T the chemical potentials. V} ^ 
is the interaction strength in the particle-particle channel, xS Mk ' (K^ Sk ') are the multipole 
(spin-multipole) interaction strengths in the p-h channel and they can be expressed via the 
Landau parameters as: 



_(M,fc) 



\ 



K 



(S,k) 



JS,k) 



K 



2r| 



/ 



/ F (r k ) ^ 

G (r k ) 
G' (r k ) 



(11) 



The monopole pair creation, the multipole and spin-multipole operators entering the 
normal products in Eq.([R]) are defined as follows: 



= E r (-i) 



j ~ m a + a + 



(12) 
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< )+ (r) = A- 1 £ T (-l) j+m (jmj' - m | A/^f (r)af m a jW , (13) 

mm 

(r) = A- 1 £ ' (-l) i+m OV - m | \l*)f$f*\T)a} m a fm ,, (14) 

j j mm 

where /./^ are the single particle radial matrix elements of the multipole and spin- 
multipole operators: 

fiS ^ u h{n)u j2 {r k )i\h\\Yx\\h) (15) 

fjih ] =u jl {rk)u j2 {r k )i L (i 1 \\T XL \\] 2 ) (16) 

One can see that the hamiltonian flTU| ) has the same form as the QPM hamiltonian with TV 
separable terms, but in contrast to the QPM all parameters of this hamiltonian are expressed 
through parameters of the Skyrme forces. 

In what follows we work in the quasiparticle representation defined by the canonical 
Bogoliubov transformation: 

a+ m = Uj a+ m + (-ly^VjOLj-rn. (17) 

The hamiltonian ( |10|) can be represented in terms of bifermion quasiparticle operators and 
their conjugates [|J: 

B(jj';\fi) = -£(-iy +m '(jmj'm | Xf^a^a^, (18) 

mm 

A + (jj';Xfi) = ^(jmj'm I x ^) a tm a } m '- ( 19 ) 

mm 

We introduce the phonon creation operators 

QU = \T. (X$ A + (jf, A/i) - (-l) x ^Yp A(jj'; X - ?)) . (20) 

33 

where the index A denotes total angular momentum and /i is its z-projection in the laboratory 
system. One assumes that the ground state is the QRPA phonon vacuum | 0), 
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i-e. Qx^i | 0) =0. We define the excited states for this approximation by Qxui I 0). For the 
QRPA the following relation is valid: 



1 



<oi [Q^otw) io) =m^ew^ -ypyp 



(21) 



The quasiparticle energies (ej), the chemical potentials (A T ), the energy gap and the co- 
efficients u,v of the Bogoliubov transformations (|T^ ) are determined from the BCS equations 
with the single-particle spectrum that is calculated within the HF method with the effective 
Skyrme interaction. Making use of the linearized equation-of-motion approach M: 



(0| SQ Xfli , H,Q+ \0)=uJO\ 6Q XfMi ,Q+ | 0) 



(22) 



with the condition: 



(o| [Qx^QtJ |o> = 



a ' 



one can derive the QRPA equations 



A B 
B -A 



\ 






fx) 






= w 













(23) 



(24) 



In QRPA problems there appear two types of interaction matrix elements, the 
matrix related to forward-going graphs and the Bfaj'^foj')^ matrix related 
to backward-going graphs. For our case we get the following expressions: 



A W 

A (iiii)r(i2i 2 ) 9 



A W 



A" 2 1 + 5 



N 



(25) 



B 



(A) 

{h>j' 1 )r{j2>j' 2 )qT 



-A 2 (l + 5 i2j < 



x 



iV 



k=l 



m 2 J no 2 



.(S,fc)\ (-) /■(AAfc) 



(-) /(AAA:). 



(26) 
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where Ejf = Ej + Ej> and uj:) = UjV^ ± VjUy. 

One can find a prescription how to solve this system and to find the eigen-energies and 



phonon amplitudes in Appendix A (see also [0). The matrix dimensions never exceed 
AN x 4A independently of the configuration space size. The derived equations have the 



same form as the QRPA equations in the QPM [0,20], but the single-particle spectrum and 



parameters of the p-h residual interaction are calculated making use of the Skyrme forces. 



III. DETAILS OF CALCULATIONS 



In this work we use generally the standard parametrization SIII |21] of the Skyrme force. 
Some examples of calculations with other parameter sets are presented in the next section. 
Spherical symmetry is assumed for the HF ground states. The pairing constants V® are fixed 
to reproduce the odd-even mass difference of neighboring nuclei. As a result constant pairing 
gaps have values that are very close to A = 12. OA -1 / 2 besides a case of semimagic nuclei. 
It is well known [12|,l|] that the constant gap approximation leads to an overestimating of 
occupation probabilities for subshells that are far from the Fermi level and it is necessary to 
introduce a cut-off in the single-particle space. Above this cut-off subshells don't participate 
in the pairing effect. In our calculations we choose the BCS subspace to include all subshells 
lying below 5 MeV. In order to perform RPA calculations, the single-particle continuum is 
discretized [2^] by diagonalizing the HF hamiltonian on a basis of twelve harmonic oscillator 



shells and cutting off the single-particle spectra at the energy of 190 MeV. This is sufficient 
to exhaust practically all the energy-weighted sum rule. As it was shown in our previous 
calculations [0 we have adopted the value A=24 for the finite rank approximation for 
the dipole and quadrupole excitations in Ar isotopes. Increasing the mass number and 
the multipolarity of excitations demands an increase of the rank to keep the calculations 
accurate. Our investigations enable us to conclude that A=45 is enough for multipolarities 
A < 3 in nuclei with A < 208. Increasing N, for example, up to A=60 in 208 Pb changes 
results for energies and transition probabilities not more than by 1%, so all calculations in 
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what follows have been done with iV=45. Our calculations show that, for the normal parity 
states one can neglect the spin-multipole interactions as a rule and this reduces by a factor 
2 the total matrix dimension. For example, for the octupole excitations in 206 Pb we need to 
invert a matrix having a dimension 2N=90 instead of diagonalizing a 1376 x 1376 matrix as 
it would be without the finite rank approximation. For light nuclei the reduction of matrix 
dimensions due to the finite rank approximation is 3 or 4. So, for heavy nuclei our approach 
gives a large gain in comparison with an exact diagonalization. 

The Landau parameters Fq, Go, F' q , G' q expressed in terms of the Skyrme force parameters 
|i~8| depend on hp. As it is pointed out in our previous work [17| one needs to adopt 
some effective value for kp to give an accurate representation of the original p-h Skyrme 
interaction. To fix the effective values of hp for the Landau parameters we use the self- 



consistency relation ||23|| . From the following set of equations: 



one can get the average field potential corresponding to such a choice of the residual inter- 
action: 



Uf f = p (to (l + ^) + ^ (h (2 + Xl ) +t 2 (2 + * 2 ))) - 



(29) 



Pr [t (- + x ) - -f (h (1 + 2x 2 ) + 2xi))j + 

-U (p a+1 (2 + a) (2 + x 3 ) - (1 + 2x 3 ) x 



24 

(2p a p T + ap a ~ 1 (pl + p 



This potential can be compared to (to* /m)U HF (r) which is the leader term of a local 
equivalent potential in Skyrme-HF approach. It is possible to evaluate the effective value k,F 
for every nucleus. One can show that this value is larger than the nuclear matter value in 
order to compensate for the effects of the neglected terms F± and G\. To calculate the dipole 

9 



strength distributions we choose kp so that the spurious isoscalar dipole mode appears at 
zero excitation energy. The strongest renormalization of the kp values in comparison with 
the nuclear matter value takes place in light nuclei. For 208 Pb the effective value kp becomes 
rather close to the nuclear matter one. 



IV. RESULTS OF CALCULATIONS 

As a first example we examine the multipole strength distributions in 36 Ar. The calcu- 
lated strength distributions are displayed in Fig.l. For the giant dipole resonance (GDR) 
and giant quadrupole resonance (GQR) QRPA gives results that are very similar to our 
previous calculations with the particle-hole RPA [[HJ] because the influence of pairing on 
the giant resonance properties is weak. This is not the case for the first 2 + and 3~ states 
that will be discussed later. For the GDR energy centroid we get E c = 19.9 MeV and this 
value is rather close to the empirical systematics 0] E c = (31. 2A -1 / 3 + 20. MeV. 



The isoscalar GQR energy centroid is equal to E c = 18.8 MeV that can be compared with 
the empirical value E c = 63 A^ 1 ^ 3 = 19.1 MeV. For the isovector GQR our calculation gives 
E c = 30.5 MeV that is about 10% less than predicted by the empirical systematics. It is 
worth to mention that experimental data for the giant resonances in light nuclei are very 
scarce. 

Results of our calculations and experimental data p5[ for the 2f state energies and 
transition probabilities B(E2) in several nuclei are shown in Table 1. One can see that there 
is a satisfactory agreement with experimental data. Results of our calculations for O and 
Ar isotopes are close to those of QRPA with Skyrme forces [0,0] and all calculations fail 



to reproduce the B(E2) value in 18 0. Making use of the SGII interaction [[Tj| improves the 



description for the O isotopes and gives practically the same results for the Ar isotopes, 
but for Sn and Pb isotopes the results become much worse. Calculations with the SkI4 
force [f2^] don't change the above conclusions. The evolution of the B(E2)-values in the Ar 
isotopes demonstrates clearly the pairing effects. The experimental and calculated B(E2)- 
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values in 38 Ar are three times less than those in 36 > 40 Ar. The neutron shell closure leads to 
the vanishing of the neutron pairing and a reduction of the proton gap. As a result there 
is a remarkable reduction of the E2 transition probability in 38 Ar. Some overestimate of 
the energies indicates that there is room for two-phonon effects. Indeed, it was found in 
calculations performed within the QPM for 208 Pb ]28|] that the two-phonon configurations 
can shift down the 2f energy by more than 1 MeV. The B(E2)-value reduction is about 10% 
in this case. The study of the influence of two-phonon configurations on properties of the 
low-lying states within our approach is in progress now. 

Results of our calculations for the 3j~ energies and the transition probabilities B(E3) are 



compared with experimental data p9| in Table 2. Generally there is a better agreement 
between theory and experiment than in the case of quadrupole excitations. Other choices 
of the Skyrme forces do not improve the agreement obtained with SHI. 

An additional information about the structure of the first 2 + , 3~ states can be extracted 
by looking at the ratio of the multipole transition matrix elements M n /M p that depend 
on the relative contributions of the proton and neutron configurations. In the framework 
of the collective model for isoscalar excitations this ratio is equal to M n /M p = N/Z and 
any deviation from this value can indicate an isovector character of the state. The M n /M p 
ratio can be determined experimentally by using different external probes fl30" H32jl . Recently 
T3| , |26[| , QRPA calculations of the M n /M p ratios for the 2f states in some O and Ar isotopes 



have been done. The predicted results are in good agreement with experimental data [26 



Our calculated values of the M n /M p ratios for the 2f and 3 1 states are shown in Tables 3 



and 4, respectively. Our results support the conclusions of Refs. [|T3|j26|] about the isovector 
character of the 2f states in 18 > 20 O and 38 Ar. As one can see from Table 4 our calculations 
predict that the M n /M p ratios for the 3j~ states are rather close to N/Z, thus indicating 
their isoscalar character. This conclusion remains valid for the SGII and SkI4 parameter 
sets. 
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V. CONCLUSION 



A finite rank separable approximation for the particle-hole RPA calculations with Skyrme 
interactions that was proposed in our previous work is extended to take into account the 
pairing correlations. The QRPA equations are derived for this case. These equations are 
used to study the evolution of quadrupole and octupole excitations in nuclei away from 
stability. It is shown that the suggested approach enables one to reduce remarkably the 
dimensions of the matrices that must be inverted to perform structure calculations in very 
large configuration spaces. 

As an illustration of the method we have used the finite rank p-h interaction derived from 
the Skyrme force SIII to calculate the energies and transition probabilities of the 2f and 3j~ 
states in some O, Ar, Sn and Pb isotopes. The values calculated within our approach are 
very close to those that were calculated in QRPA with the full Skyrme interactions. They 
are generally in a reasonable agreement with experimental data. A further development 
will be to take into account the coupling between the one- and two-phonon terms and such 
investigations are in progress now. 
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APPENDIX A: 



For the sake of completeness we show how the finite rank form of the residual forces 
(PH|) can simplify the solution of the RPA equations (pi)). In the 4A-dimensional space we 



introduce a vector 



Ps (r) 



by its components: 
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' D k p (r) ^ 



-r) 



,/3 = {M,S} 



(Al) 



where 



33 



DT (r) 



J 33 33 V 33 33 J 



33 



The index k run over the iV-dimensional space (k=l,2,...,N). Solving the system of equations 
one can get the following expressions for the phonon amplitudes: 



X% (r) = 7 r- f -=L= («« / (A , fc) + «H / (AAfc) ^ (r)) 



(A2) 



where 



y. A / (r) 

33 v 7 



E 



11 11 11 ^ 11 \ / 



+ u Xi J fc= i 



Afci V 33 33 33 33 



Xki 



2(2A + 1) : 



{D% k (r) (4 M ' fe) + 4 M ' fc) ) + ^ fc (-r) (4 M ' fe) " 4 M ' fe) )) 



2 ' 



(A3) 





(4 s ' fc) 


+ 4 s ' fe) ) 


+ (-r) 


(4^ 




^ (r) ( 


' (M,k) 




+ (-r) 


(4^ 


- 4 M ' fc) ) 



Using Eqs. (|Al|) and Eqs.([A2|),(|A3D the RPA equations fl2"ip can be reduced to the following 



system of equations: 



( Mmm (r) - 1 Mms (r) W ^ 



V 



M SM (t) A^5s ( r ) - 1 

where Ai is the 2A r x 2N matrix 



V s (r) 



(r) 



^ (Kp') + .f' fc ' ) )Tg;(r) (4' 



09',*') Jf3',k'), Tkk i / \ 



\ 



(4' 
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(A4) 



(A5) 



1 < k, k' < N. 



The matrix elements T kk ' have the following form: 

* 33 33 \ 33 J 



T k M k M (r) = T,' 



e 33 



33 



Tg 3 W 



A\k)A\\k>) (+) (-) 

/ . ./ / . ./ U . .i U . .i U!\i 
ST* JJ 33 33 33 



33 



A 2 (e 2 , 
V jj 
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Xi 



T k s k M (r) = IS5 (r) , 



« (r) = £ 



33 



AXXk) r(Wk') 

33 33 




2 

Sj' 


A 2 


V JJ 







Thus, the RPA eigenvalues u\i are the roots of the secular equation: 



det 



( M M M (t) - 1 M MS (r) ^ 



V 



0. 



/ 



(A6) 



The phonon amplitudes corresponding to the RPA eigenvalue u>\i are determined by 
Eqs.(A2), (|A3|) taking into account the normalization condition (f23|). 
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TABLES 



TABLE I. Energies and B(E2)-values for up-transitions to the first 2+ states 



Nucleus Energy B(E2|) 

(MeV) (e 2 fm 4 ) 





Exp. 


Theory 


Exp. 


Theory 


v 7 


1 Q8 


^. Id 


4^-1-9 




20 o 


1.67 


4.17 


28±2 


20 


36 Ar 


1.97 


1.91 


300±30 


310 


38 Ar 


2.17 


2.51 


130±10 


110 


40 Ar 


1.46 


2.17 


330±40 


290 


112 Sn 


1.26 


1.49 


2400±140 


2600 


114 Sn 


1.30 


1.51 


2400±500 


2100 


206p b 


0.80 


0.96 


1000±20 


1700 


208p b 


4.09 


5.36 


3000±300 


2000 
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TABLE II. Energies and B(E3)-values for up-transitions to the first 3 states 



Nucleus 



Energy 
(MeV) 



B(E3|) 
(e 2 fm 6 ) 





Exp. 


Theory 


Exp. 


Theory 


18Q 


5.10 


6.15 


1250 ±50 


1200 


20 o 


5.61 


7.28 


530 ±180 


1000 


3fi A 

d0 Ar 


4.18 


4.26 


11100 ±1100 


16000 


38 Ar 


3.81 


3.21 


10000 ±3000 


15000 


40 Ar 


3.68 


4.85 


8700 ±1000 


12000 


112 Sn 


2.36 


2.73 


87000 ±12000 


97000 


114 Sn 


2.28 


2.31 


100000 ±12000 


97000 


206p b 


2.65 


2.92 


650000 ±40000 


750000 


208p] 3 


2.62 


2.66 


611000 ±9000 


860000 




TABLE III. 


(M n /M p )/(N/Z) ratios for the first 2+ states 




Nucleus 


18Q 


20Q 


36 Ar 38 Ar 


40 Ar 


Theory 


2.1 


2.1 


0.9 0.5 


0.9 


Exp. 


0.88±0.19 a 


2.17±0.53 a 


1.41±0.50 b 


0.68±0.21 b 




1.62±0.02 c 


1.9±0.3 c 







a Ref. @, b Ref. @, c Ref. |2 
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TABLE IV. (M n /M p )/(N/Z) ratios for the first 3" states 



Nucleus 


18Q 


20Q 


36 Ar 


38 Ar 


40 Ar 


Theory 


0.8 


0.8 


0.9 


1.0 


0.9 
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FIGURES 



FIG. 1. The multipole strength distributions in 36 Ar 
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